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The surface of a 3D topological insulator is conducting and the topologically nontrivial nature of the surface 
states is observed in experiments. It is the aim of this paper to review and analyze experimental observations 
with respect to the magnetotransport in Bi-based 3D topological insulators, as well as the superconducting trans- 
port properties of hybrid structures consisting of superconductors and these topological insulators. The helical 
spin-momentum coupling of the surface state electrons becomes visible in quantum corrections to the conduc- 
tivity and magnetoresistance oscillations. An analysis will be provided of the reported magnetoresistance, also 
in the presence of bulk conductivity shunts. Special attention is given to the large and linear magnetoresistance. 
Superconductivity can be induced in topological superconductors by means of the proximity effect. The in- 
duced supercurrents, Josephson effects and current-phase relations will be reviewed. These materials hold great 
potential in the field of spintronics and the route towards Majorana devices. 



INTRODUCTION 

A three-dimensional (3D) topological insulator (TI) is a 
semiconductor with two-dimensional (2D) surface states that 
have an energy dispersion across the entire bulk bandgap. 
These surface states arise from band inversion, due to strong 
spin-orbit interactions. The band inversion makes the surface 
states topologically nontrivial, meaning that these states are 
protected (cannot be removed) and that the states have a he- 
lical Dirac-type dispersion, where spin is tightly coupled to 
momentum. Several compounds have been theoretically pre- 
dicted to be 3D TIs |1|. Strong experimental evidence for 
the presence of topological surface states exists for the class 
of Bi-based materials: the alloy Bii-jjSb^; fT\, and the com- 
pounds Bi2Se3 and Bi2Te3 [31, see Fig. [T]for the crystal struc- 
ture. The class of Bi-based topological insulator has been ex- 
panded with the fabrication of TlBiSea f4^-Vi, TlBiTca 0, 
BiSbPbTe |7|, andPbBiTe4 [8|. The bulk of these compounds 
is not necessarily insulating, due to the presence of defects 
and impurities. This bulk conductivity might shunt and mask 
topological transport properties. A successful trend in the re- 
search on topological insulators has been to reduce the num- 
ber of defects and impurities by atom substitutions such as 
Bi2Te2Se|9-ll I and Bi2-j;Sba;Te3_j^Sej, |12|. which have 
higher bulk resistivities. Simultaneously, progress has been 
made in disentangling topological surface state transport prop- 
erties from bulk transport contributions. Here, we review and 
analyze the experimentally observed magnetotransport prop- 
erties of the Bi based 3D TIs as well as the induced super- 
conductivity in these compounds. Studying transport proper- 
ties in an externally applied large magnetic field is a powerful 
probe of the topological nature of the surface states. A Berry 
phase shift appears in magnetoresistance oscillations and in 
quantum corrections to the conductivity, leading for example 
to weak anti-localization effects. We will review the ongoing 
debate on the manifestation of this phase shift and show how 
the phase shift depends on the actual parameters of the sample. 
We will discuss the important role of bulk conductivity contri- 
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FIG. 1: The crystal structure of Bi2Te3 as a representative example 
of Bi-based 3D topological insulators. From jS). 



butions. Intriguingly, two-dimensional transport features not 
always originate from the surface in Bi-based topological in- 
sulators. This finding is also important for understanding the 
presence of a large and linear magnetoresistance in most Bi- 
based 3D TIs. 

The topologically non-trivial surface state is an attractive 
medium to study unconventional superconductivity. When 
superconductivity exists in a system with a Dirac-type dis- 
persion and helical spin-momentum locking, exotic features 
could be realized such as p-wave order parameter symme- 
try |fT3]| and Majorana type bound states |14|. Josephson 
hybrid structures have been designed theoretically already to 
mimic non-abelian particle statistics 1 14 |. Parallel to the trend 
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of doping and alloying TIs to render them superconducting 
ifTSi [T6ll runs the successful approach of inducing supercon- 
ductivity from a standard superconductor into a TI by means 
of the proximity effect. Here, we will review the current sta- 
tus of the experimental efforts in the latter direction. Evidence 
for supercurrents through topological surface states will be an- 
alyzed, also in the presence of bulk conducting channels. The 
prospects for future SQUID based experiments to detect topo- 
logical effects in the superconductivity will be discussed. 



WEAK ANTILOCALIZATION 

Weak antilocalization (WAL) is commonly observed in 
magnetoresistance measurements on thin films or flakes of Bi- 
based 3D TIs, see Fig. |2|a). WAL is observed for a wide 
range of thicknesses. Kim et al.fVf] report magnetotransport 
on Bi2Te3 films from 3 nm to 170 prn thick, where they can 
still see the WAL effect. When caused by topological surface 
state transport, WAL may more properly be named topologi- 
cal delocalization lfTSll , because the antilocalization is caused 
by the tt -Berry phase (see section 3) picked up along two 
time-reversed self-intersecting paths. Effectively, in TI sur- 
face states all the WL orbits contribute as WAL orbits. WL 
will not appear in TIs, as long as time reversal symmetry is 
not broken. 

The Hikami-Larkin-Nagaoka (HLN) equation |fT9l de- 
scribes the quantum correction to 2D conductivity due to 
(anti-)localization effects as 



^(^xx{B) = (Jxx{B) - crj;^(0) 



~B 



In 



(1) 



where ^' is the digamma function, = j^^jj, and is 

4' 

the dephasing length. Strong spin-orbit interactions (or topo- 
logical effects) give WAL with a = —0.5 while weak spin- 
orbit scattering typically gives weak localization (WL) effects 
with a = 1. A logarithmic temperature dependence is ex- 
pected, since L^, generally scales inversely with temperature. 
For WL the conductivity decreases with decreasing tempera- 
ture and for WAL the conductivity increases with decreasing 
temperature. An applied magnetic field will reduce the effect 
of the quantum corrections, due to a phase being picked up in 
the magnetic field, destroying the interference effects. Equa- 
tion ([T]) is an approximation of an expression that takes all the 
scattering lengths explicitly into account IJ9i . 

For topological insulator devices, the value of a in Eq. ([T]i 
is usually around —0.5, but it is found to depend on sample 
thickness, the presence of bulk conductivity, gating and the 
substrate, as will be elaborated below. The trivial value of 
a = is only found for higher temperatures or for very thin 
samples in a parallel applied field 1,20] . Depending on the 
number of conduction channels, a can even become —1 when 
both the top and the bottom surface of a crystal or film con- 
tribute independently [;21h23 J . The 2D nature of the WAL is 



confirmed by angle dependent measurements. For a 2D sys- 
tem, the WAL effect only depends on the perpendicular com- 
ponent of the field 1201 l24l . For higher fields extra spin effects 
may arise due to the Zeeman splitting caused by the parallel 
field I 



Bulk vs surface 

Even when the 2D nature of the WAL is confirmed by angle 
dependent measurements, the WAL can still be caused by bulk 
effects. Quite a number of reports U2\ \M EOJEI |27H33]| 
deduce a dephasing length from the WAL fit using the HLN- 
equation that is larger than the film thickness, making the bulk 
effectively more 2D. Also, it is an intriguing question why in 
most cases the observed WAL can only be explained by just 
one of the two surfaces (i.e. a — —0.5). One option would be 
that the whole topological insulator can be considered as one 
2D system with spin-orbit coupling, as a mix of bulk channels 
and surface states Il20lj22l 1311 132 1 . The question then is how 
these contributions can be disentangled. An answer might be 
found in the hints of bulk contribution found in some works, 
such as a broader WAL from the bulk [24], WL channels 134}- 
l36l or by specifically trying to influence the surface states, 
for example by gating 12114231 133 1 or by introducing magnetic 
particles |20|. A theoretical description of the bulk-surface 
coupling is given by Refs. |35| and |36|. 

He et al. 1201 show that the WAL in Bi2Te3 thin films has 
contributions from both the surface and the bulk, by analyz- 
ing angle dependent measurements. The surface contribution 
to the WAL can be suppressed by depositing magnetic impu- 
rities (Fe) on the surface. For very thin films they observe a 
crossover to the regime of a field dependence. The ques- 
tion does arise what happens to the bottom interface; the film 
is smaller than the phase coherence length and the bottom 
interface should not be affected by the magnetic impurities 
on the top interface. The authors note that the contributions 
from the bottom surface seem negligible and that this may be 
caused by defects present at the Tl-substrate interface. Check- 
elsky et al. [23 1 find contributions from two interfaces us- 
ing the HLN equation to fit the WAL measured on exfoliated 
flakes of Bi2Se3, while still having a phase coherence length 
larger than the thickness of the flakes. This gives merit to the 
claim that for grown thin films, the bottom interface has too 
much disorder to show a WAL effect |20|. 

Steinberg et al. Ii22ll observe a tunability of the number 
of channels leading to WAL in top gate tuned 20 nm Bi2Se3 
thin films, both by temperature and applied electric field. The 
authors argue that the different channels can be viewed as in- 
dependent as long as carriers in one channel lose coherence 
before scattering to the other channel, i.e. tsb > t^, where 
TsB is the surface-to-bulk scattering time. This can be the 
case when a depletion layer is formed between the surface and 
the bulk by applying a gate voltage l2Bl23l [3Tll33ll . Stein- 
berg et al. 1221 observe that the HLN fit also works in the 
crossover regime between a = —1 and —0.5, indicating that. 





0.75 - 



1.4 



T — 


-1 — 1 — 1 — ^ 
/•• • 


;> 


,0.,.;;-» . 






1 


I.I. 



2 4 
B (T) 



0.3 1 3 

Temperature (K) 



10 



x=0.13 




x=0 (xO.5) 



-10 12 



FIG. 2: Weak antilocalization in topological insulator thin films and crystals, (a) thickness dependence of WAL for thin Bi^Ses films by Liu 
et al.\ 18], the red line is a fit with the HLN equationi 19|, Eq. l[TJ. (b) Temperature dependence in magnetic field of an 80 nm Cua:Bi2-a;Se3 
film (x = 0.02 - 0.03) 1281 . The logarithmic temperature dependence is largely field independent. The inset shows / and Tq determined from 
a fit to Eq. where / = (l - |F) and To = T^e. The WAL contribution is suppressed above 2 T. (c) WAL to WL crossover observed in 
Bi2_a,Cri,Se3 films as a function of Cr doping 1401 . Red lines are fits using the HLN equation. 



while there is no theoretical model, a can be used as a phe- 
nomenological measure for the channel separation (i.e. the 
ratio tsb/t^)- As the temperature is increased the channel 
separation is seen to increase, which is explained by the strong 
temperature dependence of versus the relatively tempera- 
ture independent tsb^ which is suggested to be governed by 
impurity scattering 1371 . A similar empirical model is sug- 
gested by Kim et al. ifTTll . In their measurements on up to 
100 unit cells of Bi2Se3, the mobility depends linearly on the 
film thickness, indicating that should also depend linearly 
on film thickness if the WAL would be a purely bulk effect. 
However, in the thin film limit the scaling of deviates from 
a linear scaling, indicating it cannot be explained by only a 
bulk contribution. The authors note that the scaling of 
with thickness t can be used as a figure of merit for TIs. For 
(X t^, s should be zero for a perfect TI with insulating 
bulk, one for a trivial strong spin-orbit coupling material or 
somewhere in between for a TI with both surface and bulk 
contribution. 



Localization to antilocalization crossover 

A combination of WL and WAL is observed at low tem- 
peratures in ferrocene doped Bi2Se3 nanoribbons ll24l . This 
can be explained by the opening of a band gap in the surface 
states induced by the magnetic impurities |38 1 or by a cancel- 
ing of the bulk spin-orbit coupling by the magnetic impurities 
E4\ . A more detailed study with controlled magnetic doping 
is needed to disentangle these two effects. A coexistence of 
WL and WAL for nonmagnetic dopants has not been observed 
Il20ll24l|39 !l, which is an indication that (a significant portion 
of) the WAL is caused by the topologically protected surface 



states. 

Liu et al. I;40J observe a crossover as a function of Cr dop- 
ing in Bi2Se3 films of 3 unit cell quintuple layers (QL) grown 
by molecular beam epitaxy on sapphire, see Fig. |2|c). While 
this thickness is already in the regime where scattering be- 
tween the two surfaces opens up a gap in the surface states 
|41 1, the authors note that the observed effects are still nomi- 
nally the same as in 6 QL films, but without a quadratic bulk 
contribution. The crossover looks similar to what has been 
predicted theoretically 1 38^, driven by the opening of a gap due 
to the magnetic impurities Il42ll43]| . However, the authors note 
that from their ARPES measurements, it looks like the surface 
states completely disappear for large magnetic doping, leav- 
ing the system in a dilute magnetic semiconductor state. They 
find a — 0.14 for the HLN fit for the highly doped films, 
which is indeed closer to the regime of strong magnetic scat- 
tering than the WL regime |,19J. The authors observe a further 
crossover as a function of temperature, where for the highly 
doped samples there is a coexistence of WL and WAL, each 
with a different critical field scale 1,401 . This coexistence could 
be another explanation for the low value of a; in the crossover 
regime the limiting cases used to obtain Eq. ([T]i do not hold 
and the full expression should be used L19..44J . 



Electron-electron interactions 

Electron-electron interactions (EEI) can lower the conduc- 
tivity via the Aronov-Altshuler effect Il45ll46l . The Coulomb 
interaction between different electrons decreases the effective 
density of states near the Fermi energy, which results in a cor- 
rection to the conductivity that is logarithmic in temperature. 
In magnetic fields the conductivity will further decrease as a 
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result of Zeeman splittingf46'|. 

The correction to the conductivity by EEI |46 | is given by 

where F is the Coulomb screening factor and T^e is the char- 
acteristic temperature for EEI effects. 

EEI is used by some authors ITSl l28l |47| to explain a de- 
crease in conductance that depends logarithmically on tem- 
perature, where an increase is expected for just WAL. This 
decrease is not always observed: in some cases there is only 
a metallic temperature dependence 1211 [32l . EEI have already 
been used to explain localization effects in other strong spin- 
orbit materials |48 , 49 1. The EEI contribution to conductivity 
is extracted from the conductance in magnetic fields. 

Checkelsky et al. Il47l deduce both orbital and spin contri- 
butions to the magnetotransport from the angle dependence of 
universal conductance fluctuations. The low field effect is fit 
by a combination of WAL ||T9| (the orbital part) and EEI ||46| 
(the spin part). They find a logarithmic field dependence, but 
the fitting parameters are orders of magnitude larger than ex- 
pected for a 2D system. The authors reason that extra 2D bulk 
states with similar spin-orbit properties as the surface states 
may be able to explain this discrepancy [50]. 

Takagaki et al. |28 1 determine the prefactor and T^e in Eq. 
Q by looking at the temperature dependence in magnetic 
fields, to suppress the contribution of WAL, see Fig. |2|b). 
The value they find for the screening factor is too large to only 
come from the surface, which suggests that the bulk states also 
have to be taken into account. 

Liu et al. |[18 1 observe a combination of WAL and EEI in 
ultrathin Bi2Se3 films of one to six QL grown on sapphire. 
The logarithmic temperature dependence is determined to be 
caused by EEI by studying the MR in parallel fields. EEI 
is largely field independent, but for large Zeeman splitting it 
will be suppressed with a logarithmic field dependence. The 
g-factor is not very large in TIs ISTl . but Liu et al. indeed 
observe a logarithmic field dependence for fields above 4 T. 
They argue that the EEI becomes larger for thinner films due 
to increased disorder and reduced dimensionality, leading to 
poorer screening and therefore stronger Coulomb interaction 

SHUBNIKOV-DE-HAAS OSCILLATIONS AND BERRY 
PHASE 

The mobility fi of the Bi-based topological insulators has 
frequently been reported to be high enough to fulfill the con- 
dition fiB ^ 1, where Shubnikov-de-Haas oscillations can be 
observed in the resistance at high magnetic fields 1 10 , 32, 51- - 
l63l . The oscillations in the longitudinal surface resistance ^ 
are given 1641 by 

Apf, cx cos (2^^+71 + 0s) , (3) 



where (pB is a Berry phase, as will be discussed below, and the 
frequency F is determined from the extremal cross-section of 
the Fermi surface, 2ttF = hS{Ep)/e. 

If the electron orbit is confined to a 2D plane as for a two- 
dimensional electron gas, then B is the perpendicular mag- 
netic field component. However, a cosine magnetic field an- 
gle dependence of the oscillations alone is not enough to con- 
clude that the topological surface states are providing the os- 
cillations. For example, it has been shown for Bi2Se3 that the 
bulk may act as the sum of many parallel 2D electron systems 
|61 1. And a surface state might also be a non-topological triv- 
ial band that crosses the Fermi energy due to band bending at 
the surface, as observed for Bi1.5Sbo.5Te1.7Se1. 3 1581 . 

Rather, the topological nature can be concluded from the 
phase of the Shubnikov-de-Haas oscillations. From graphene 
it is known how a Dirac cone gives rise to a half filled Landau 
level at zero energy (here defined as the energy of the Dirac 
point) and how this Landau level gives an additional Berry 
phase shift of tt to the quantum Hall effect or Shubnikov-de- 
Haas oscillations ||55l IS6l . In general, the Berry phase can 
be obtained by integrating the Berry connection fl over an 
electron orbit 




where the Berry connection follows from a specific band 
structure as 

j <(r)VfcMk(r)dr, (5) 

where Uk(r) are the periodic parts of the Bloch wave func- 
tions. It turns out that Eq. Q gives (t)B = tt at the surface 
of the Bi-based topological insulators irrespective of the ex- 
act shape of the Dirac cone l67l . i.e. deviations of the linear 
dispersion in a Dirac cone do not lead to different values for 
the Berry phase, in contrast to previous expectations l68ll . The 
crucial ingredient to have a tt Berry phase is that the order of 
conduction band and valence band are inverted due to strong 
spin-orbit coupling. 

For the trivial (t>B = case, Eq. ([3]) shows that the 
magnetoresistance oscillations have a minimum in the limit 
1 /B — > 0. A plot of the 1/B positions as function of the n* 
number of minima then extrapolates to 0. When 0^ = tt, 
the 1/B positions of the minima extrapolate to —1/2. The 
extrapolation becomes very accurate either when the applied 
magnetic fields are so large that the lowest Landau levels are 
resolved, or when an additional gate voltage is used to sweep 
the Fermi energy [55 1. 

Inverting the resistivity tensor gives the longitudinal con- 
ductivity as axx = Pxx/ {Pxx + Pxy)- In ^ one-band free 
electron model, the oscillation condition /iS ^ 1 coincides 
with pxy > Pxx since p^x = {enp,y^ and p^y = {eny^B, 
where n is the carrier density. When the oscillations in p^x 
are small {/S.pxx ^ Pxx) and when oscillations in p^y can be 
neglected, the fulfilled condition p^y > Pxx automatically en- 
sures that A(Txx oc Apxx, i-e. the oscillations in conductivity 
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and resistivity are in phase. Both a^x and pa:x have minima 
at integer Landau level index n. Magnetic field measurements 
on semiconductors are usually in this regime. 
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FIG. 3: Landau level index n as function of l/B for Bi2Tea and 
Bi2Te2Se. (a) For thin flakes of Bi2Te3 the Shubnikov de Haas oscil- 
lations are shunted by a bulk conductivity with a moderate mobility. 
In this case, minima positions (black squares) of p^^a; coincide with 
integer n and extrapolation to l/B = provides a Berry phase of tt, 
indicating the topological nature of the surface states. Data is taken 
from Veldhorst etal. |60|. (b) For a large crystal of Bi2Te2Se with a 
bulk conductivity shunt with a low bulk mobility, the integer values 
for n appear for the maxima positions (filled symbols) of p^^^. The 
inset shows the Landau level filling at 42 T (arrow). From 1511 . 

However, in the presence of a buUc shunt with a low mobil- 
ity, the oscillations in p^-j. and cr^x can be out of phase. It was 
shown for a Bi2Te2Se crystal showing bulk conductivity with 
a low mobility that the surface state resistivity has maxima 
now at integer n 151]. Consequently, the Berry phase needs to 
be inferred from a plot of the XjB maxima positions at integer 
n, see Fig. [3|b), where the topologically non-trivial 4>b — 
was observed. However, this situation L5 1 J is not generic for 
any strong bulk conductivity shunt. In a parallel summation 
of surface and bulk conductivity, it can be shown that a^x and 
Pxx are in phase when the bulk mobility is not too low, i.e. as 



long as tr^ 



> at 



For example, Cao et al. [61] have shown 



that the bulk in highly doped Bi2Se3 can act as the sum of 
many parallel layers with high mobility and a trivial Berry 
phase with minima at integer n. Lower bulk densities will en- 
hance the surface state conduction. In Fig. [Sj^a) data is shown 
for a 100 nm thick exfoliated flake of Bi2Te3. A multiband 



analysis (assuming two surfaces and one bulk channel) of the 
Hall data and Shubnikov de Haas oscillations of Ref. |l60| 
provides a bulk carrier density of 1.9 x 10^^ cm^'^ and a bulk 
mobility of the order of 0.1 m^V~^s~^, so that a^y > cr^j, is 
fullfilled above about 10 T. The extrapolation of the minima in 
Pxx at integer n consistentently extrapolate to (ps = tt again. 
For every sample with bulk conductivity a detailed parameter 
estimate should decide in which regime the oscillations are 
and whether maxima or minima positions should be consid- 
ered. 



LINEAR MAGNETORESISTANCE 

In addition to quantum transport phenomena such as 
Shubnikov-de Haas oscillations and weak antilocalisation, 
a strong background magnetoresistance (MR) was also ob- 
served in the 3D TIs BiaTeg l52l|69l, Bi2Se3 ED [33] Ei ED, 
and Bi2Te2Se 11721 . The reported magnetoresistance varies 
from a few percent to more than 600% [69 J . In contrast to 
the more standard small quadratic MR (e.g. due to multiband 
effects), which saturates for fiB > 1, the observed MR is 
linear with field and does not seem to saturate in high fields. 
The MR is positive, also for negative field directions, thereby 
ruling out Hall component admixtures due to misaligned elec- 
trodes. 

Linear MR is not a new phenomenon and was already ob- 
served about 100 years ago in Bi 1731 [741 . The so-called 
Kapitza linear MR as observed in metals with open Fermi 
surfaces require large Fermi areas (and therefore large car- 
rier densities) and this does not really apply to any of the Bi- 
based 3D TIs of current interest. Moreover, more quantita- 
tive treatments of open orbits seem to indicate B^^"^ or i?"'/^ 
rather than linear behavior. Abrikosov proposed a quantum 
MR model to account for linear MR 1,74,-761 when only one 
Landau magnetic band is filled. The conditions for the sole 
occupation of the n — lowest Landau level (at energy Eq) 
is that both the Fermi energy, Ep, and ksT are much smaller 
than Ei^Eq, where Ei is the energy of the next Landau level. 
This is called the extreme quantum limit and the resulting lin- 
ear MR is called quantum linear MR. 

Usually, the conditions for the extreme quantum limit are 
only fulfilled for narrow gap semiconductors or semi-metals 
with very small Fermi pockets and low effective masses. Ex- 
amples are Bi, n-type doped InSb llTTl . and PbS ifTSj . For a 
parabolic dispersion, the Landau levels are given by En = 
(71+5). A low effective mass gives a large cyclotron 
frequency, lUc — ^ and a large spacing between Landau lev- 
els. Entering the extreme quantum limit is also facilitated in 
materials with a linear dispersion. E.g., for /3-Ag2Te a lin- 
ear MR was observed for fields as low as 10 Oe at 4.5 K 
Il79l . for which the conditions of the extreme quantum limit 
do not seem to be fulfilled (fc^T ^ TluJc for B = 10 Oe 
and TO = O.OlTOe gives T < 0.1 K). However, for a ma- 
terial with linear dispersion, the Landau levels are given by 
En = vnV'^ehBn, where is the Dirac velocity. This pro- 
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vides a less stringent condition for the extreme quantum limit 
at 10 Oe (T < 4 K for w^, = 3 x 10^ m/s). Recendy, it has 
been calculated that /3-Ag2Te is indeed a material with a linear 
dispersion |80|. 

Apart from a single Landau level explaining linear MR, 
also inhomogeneities in disordered conductors can give rise 
to linear MR. This is described by the model of Parish and 
Littlewood ifSTl [82l . Regions with higher conductivity create 
deviations of the direction of the applied current and classi- 
cal resistor network models show that also in this case a Hall 
voltage admixture occurs, giving rise to a linear MR Ii8n [82l. 
The Hall signal mixes in such a way that the MR is still sym- 
metric in field. The magnitude of the magnetoresistance from 
this effect will not exceed the maximum local Hall resistance 
in the inhomogeneous system. Note, that also the classical in- 
homogeneity model was used to explain the /?-Ag2Te linear 
MR |,81J . 

The nature of the linear magnetoresistance in topological 
insulators is not clear yet. In the past few years a lot of 
research has focused on unraveling the mechanism behind 
this linear magnetoresistance and to investigate which model 
applies and whether or not linear MR is a signature of the 
topological surface states. Furthermore, unravelling the ori- 
gin of the linear magnetoresistance can have important con- 
sequences for spintronics, as the observed magnetoresistance 
is already large in amplitude and survives even up to room 
temperature. 

Bulk vs surface 

Qu et al. |52| found a magnetoresistance of 1.7% at 0.3 K 
in a magnetic field of 14 T in Bi2Te3 crystals. The observed 
dependence on perpendicular field is linear, starting already 
from a very small magnetic field of 30 G. The in-plane MR 
is also linear, but the effect is a factor of 2 smaller than for 
perpendicularly applied field. The authors attributed the MR 
to the coupling of the spin of the bulk electrons to the mag- 
netic field. It is interesting to note that the linear magnetore- 
sistance (LMR) in this case was only observed in the non- 
metallic samples, while the metallic sample showed quadratic 
MR. Furthermore, besides the LMR also Shubnikov-de Haas 
oscillations were observed for the out-of-plane field. 

The angular dependence of the LMR was investigated fTOl 
ITTI . By changing the direction of the field from out-of-plane 
to in-plane, the LMR decreased significantly and a | cos(6')| 
relation was found between the LMR and angle of the field. 
This indicates that the LMR has some 2D character. Tang et 
al. |70 | deduced the position of the Fermi energy from the 2D 
sheet carrier density from the SdH oscillations and found this 
to be consistent with surface states, althought the expected 
TT Berry phase could not be clearly resolved. He et al. lITTl 
found that the LMR disappears for ultrathin films. This could 
be related to the gap opening of the surfaces state due to the 
coupling of the two surfaces with each other 141]. 

While it is tempting to attribute the LMR to the topological 



surface states, as bulk transport generally has a 3D character, 
Cao et al. \6l \ have shown in their highly doped Bi2Se3 sam- 
ples that the observed 2D Shubnikov-de Haas oscillations can 
originate from the bulk. It is possible that the individual quin- 
tuple layers in the bulk all act as 2D parallel electron systems. 
Solely observing the | cos(6')| is thus not sufficient to attribute 
the magnetotransport effect to the surface states. 

Gao et al. investigated the MR as function of electric gate 
voltage in Bi2Se3 sheets ll33l . They found that by increasing 
the negative gate voltage, the quadratic MR changes to a LMR 
and also the effect of WAL enhances. Comparing this to the 
temperature dependence of the resistance, a peak appears in 
the metallic temperature dependence by increasing the gate. 
This indicates that two components contribute to the charge 
transport. The thermal excitation of carriers is attributed to 
an impurity band in the bulk and the metallic component to a 
state in the gap. It was also assumed that the impurity band 
gives a quadratic MR, while the gap state has a linear con- 
tribution to the MR. Fitting of the MR and the temperature 
dependence results in consistent values for the gap and bulk 
state resistance values. By fitting the WAL effect with the 
HLN equation, Eq. ([TJ, they found that a increases from -0.6 
to -1.2 with increasing gate voltage. Therefore, they attributed 
the metallic component to the topological surface state. How- 
ever, again it is difficult to distinguish whether the WAL is 
due the 2D surface states or that the whole sample acts as a 
2D system, as is noted in section 2. 
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FIG. 4: High field measurement of the magnetoresistance of a thin 
flake of Bi2Te3. The field was applied from out-of-plane (top curve) 
to in-plane (bottom curve) in steps of 10 degrees. The inset shows a 
cosine fit to the data. This dataset was published in Ref. [601 after 
substracting the linear backgrounds. 

An important observation in the discussion whether 2D 
effects originate from the bulk or the topological surface is 
shown in Fig. |4] The dependence of the MR on the applied 
magnetic field direction suggests 2D behavior (although the 
MR does not vanish for in plane fields, see also Ref. [52 1). 
However, from analyzing Hall measurements and Shubnikov 
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de Haas oscillations, it can be shown that the conductance of 
the bulk of this sample is larger than the conductance from the 
surface 1601 . Therefore, when a large 2D effect is then seen 
on the scale of the total conductace, it seems unlikely that the 
LMR in this case originates from the surface states alone. 

Quantum vs classical effect 

Attempts have been made li69]j7TI to explain the obser- 
vation of LMR by the quantum magnetoresistance model of 
Abrikosov. The disappearance of the LMR for ultrathin films 
indicates a relation to the gapless linear dispersion. The low 
temperature dependence of the observed LMR also points to- 
wards Abrikosov's model. However, as Abrikosov's model 
applies to the extreme quantum limit, it is interesting to note 
that not all LMR observations |'69"7T1 coincide with observa- 
tions of Shubnikov-de Haas oscillations. This could indicate 
that the systems are far away from the extreme quantum limit 
or that the bulk transport shunts the quantum oscillations. The 
carrier density in the topological insulator samples are rela- 
tively large and for realistic magnetic field values more than 
one Landau level is usually occupied. If Abrikosov's model 
is used to explain the LMR in these topological insulators, it 
is important to know up to which filling factor Abrikosov's 
model could still apply. 

To explain the LMR observed in Bi2Se3 nanosheets ||70l , 
Wang and Lei |83| derive a different model for the LMR us- 
ing the balance-equation approach 1(841 . They found that for 
a 2D system with linear dispersion and a non-zero g factor 
a LMR can arise. The model shows that the LMR can oc- 
cur for fields much lower than iiB — 1, it actually requires 
that the Landau levels overlap. When the mobility decreases, 
the required field for SdH oscillations increases, whereas the 
LMR gets even stronger For higher mobilities, SdH oscilla- 
tions start to appear and the LMR weakens, because in this 
regime the overlap of the Landau levels decreases. This also 
accounts for the low temperature dependence of the LMR, as 
with higher temperature the Landau levels will broaden, in- 
creasing the overlap. It is clear that the system is far away 
from the extreme quantum limit, in contrast to Abrikosov's 
model. By investigating the density dependence of the LMR, 
the model of Wang and Lei can also be distinguished from 
Abrikosov's model. The latter predicts an R^x oc N^'^ be- 
havior where s — 2, while the former predicts s = L A 
positive g factor is required for a positive LMR. The question 
arises whether Wang and Lei's model could also explain the 
strength of the LMR, as some reports indicate that the sur- 
face g factor could be quite low [51 1. This model was derived 
for Bi2Se3 nanosheets |70|, but could deviate for other topo- 
logical insulator samples. It assumes a linear dispersion, but 
ARPES measurements have shown that Bi2Te3 deviates from 
a linear dispersion. The bulk contribution was assumed to be 
negligible, although in many transport measurements the bulk 
is still prominent. It would be interesting to see whether the 
model by Wang and Lei [83 ] would still be linear when gener- 



alized beyond a linear dispersion and in the presence of bulk 
shunts. 

Alternatively, Assaf et al. Il72l used a modified HLN model 
to explain the observed LMR. Their extended HLN model 
takes the WAL effect and a classical component into ac- 
count. In a certain magnetic field range, the logarithmic WAL 
component cancels the component, leaving a linear MR 
dependence. According to this model, the quadratic com- 
ponent should be distinguishable from a truly linear MR for 
higher fields, so it would therefore be interesting to investi- 
gate whether this idea is still consistent with measurements at 
even higher fields. 

A number of reports (|33l |69l |69] |70l EH UlSl mention that 
the Parish-Littlewood model also could cause a LMR. Al- 
though Assaf et al. ||72][ indicate that their samples contain 
a lot defects, they exclude the possibility that this model ex- 
plains the LMR, as there is a disagreement between the range 
of the mobility and the range of the LMR. Others [[69] EH 
also indicate that it is unlikely that the model by Parish and 
Littlewood explains the observed LMR, as the topological in- 
sulator samples have a high quality crystal structure and no 
structural inhomogeneities are observed. Although Tang et 
al. [70[ also do not expect structural inhomogeneities, they 
indicate that electronic inhomogeneities cannot be excluded, 
which has been confirmed by the observation of potential fluc- 
tuations in scanning tunneling measurements [85 J . More re- 
search is therefore required to investigate whether the model 
by Parish and Littlewood could apply to these topological in- 
sulator samples. Expecting a relatively large size of electrodes 
or contacts in some of the reported transport experiments it 
would also be useful to study more mundane explanations 
of LMR such as a non-homogeneous current injection into 
an otherwise homogenous sample, whereby a transverse Hall 
component can easily be picked up li86l . 

TOPOLOGICAL INSULATORS WITH PROXIMITY 
INDUCED SUPERCONDUCTIVITY 

The combination of superconductivity with spin-orbit cou- 
pling opens many new exciting research directions, includ- 
ing the realization of a new emergent particle: the Majorana 
fermion. This particle, which is its own antiparticle, emerges 
from the electron-hole symmetry of superconducting conden- 
sation and the spin-momentum locking in the topological in- 
sulator lfT4l [87l . The high potential of this particle in quan- 
tum computation led to many new proposals in various ma- 
terial systems, such as topological insulator/superconductor 
structures lfT4l [88l [89 [, spin-triplet p-wave superconductors 
lfT3ll9QH93l , helical superconductors [94[, topological super- 
conductors lfT5l[T6ll95l j97l, semiconductor/metals with strong 
orbit coupling in combination with superconductors or super- 
conductors with strong spin orbit coupling [,98l4101L In the 
latter category the most well-known device is inducing su- 
perconductivity in nanowires structures with strong Zeeman 
and Rashba fields MlOOllTTl j . First signatures of the Majorana 
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Current (|jA) 

FIG. 5: Bessel-peacock color plot of the conductance after irradiation 
with microwaves (10.0 GHz). From L60J . 



shorter than the junction length while the system is in the bal- 
listic limit and concluded that the discrimination between bulk 
and surface states can be so strong that the top topological sur- 
face state dominates the superconducting transport, while the 
bulk states dominate the normal state. 

It is an intriguing question why this discrimination is so 
strong. Besides geometrical effects, strong anisotropy of the 
band structure of Bi-based topological insulators is a pos- 
sible origin. The ratio of the superconducting coherence 
length ^zKxy in the topological insulator depends on the ef- 
fective mass and the potential energy difference along the 
xy-direction and z-direction. From the expression of the 
coherence-length tensor for anisotropic superconductors as 
described in Ref. BI23I . we find for the proximity coherence 
length in the non-superconducting material 



fermion are revealed HI 021 11031 1112H114I but in order to test f m h m 

all its peculiar properties the realization of superconducting 7 — cx: — oc 

interference devices will be a prerequisite 111 15111 16L 



Supercurrents through TI structures 

Shortly after the discovery of the Bi-compound topolog- 
ical insulators, observations of supercurrents in supercon- 
ductor - topological insulator structures have been reported 
|55 60 1 171 412111 . as well as the coexistence of superconduc- 
tivity and topological surface states 1 122|. Several fabrication 
techniques have been used such as epitaxial film growth and 
mechanical exfoliation for the topological insulator structure. 
Superconducting electrodes are deposited on top resulting in 
lateral junctions. Sandwich-type junctions have been realized 
by Qu et al. IIII7I using a clever fabrication method based on 
exfoliated flakes. Using an oxidized Si substrate, Sacepe et 
al. |I55| demonstrated gate tunable supercurrents. The change 
in the IcRn product is of the order of 15% by applying a gate 
voltage between -80 V and H-50 V. Structures based on exfoli- 
ation are therefore flexible and have high potential. High qual- 
ity topological insulators with small thicknesses and smooth 
surfaces can be realized that allow the construction of multiple 
devices on a single flake |60|. 

The characteristic IcRn product observed so far is well be- 
low the characteristic voltage scale of 7rA/e, e.g. 5 mV for 
superconducting niobium electrodes. Although mechanisms 
as scattering and electron-hole decoherence due to junction 
lengths longer than ^ can lower the IcRn product, these junc- 
tions often have high interface transparencies and junction 
lengths comparable to ^. In these systems, an important mech- 
anism lowering the IcRn product is the presence of bulk con- 
duction. Normal state transport is dominated by the bulk states 
which are spread over the entire sample. As the electrodes are 
fabricated on top, proximity induced superconductivity will 
be the strongest on the top, and the bulk is effectively shunt- 
ing the junction Rn significantly, thereby reducing the IcRn 
product. Veldhorst et al. [60] obtained a bulk mean free path 



where the x, y and z axis are as defined in Fig. [T] The super- 
conductivity is induced in the surface states in the xy— plane. 
a and c are the lattice constants of the unit cell along the x 
(or y) and z direction, respectively. Using the valence band 
effective masses rrij-y = 0.18 and = 0.84 for Bi2Te3 
111 24 J in conjunction with the lattice constant differences, we 
find a ratio (^zji^xy of 0.15. Hence, the band structure itself 
causes a strong anisotropic proximity effect. An additional 
argument for strong anisotropy is the observation of the lay- 
ered quantum Hall effect in Bi2Te3 1611 . Furthermore, the 
effective parameter that determines the suppression of IcRn 
in S-N-S junctions (whether ballistic or diffusive, and includ- 
ing barriers at interfaces) is given by 7e// — (Jf)'^^^^ with 
the boundary resistance Rb and the interlayer resistance Rn 
ltl25J . From this suppression parameter it can akeady be seen 
that a low interlayer resistance with respect to the boundary re- 
sistance gives a high 7e//- The IcRn is then suppressed more 
strongly for the bulk than for the surface channel. Finally, 
topological effects might render the interfaces to the surface 
states intrinsically transparent. These intriguing options open 
a new exciting research direction and give the opportunity to 
have surface states dominating the supercurrent in topological 
insulators even in the presence of substantial bulk conductiv- 
ity. 

Successive Andreev reflections allow quasiparticles to es- 
cape from the superconducting gap in the voltage state. These 
multiple Andreev reflection (MAR) processes cause structures 
in the cuiTent voltage characteristics at voltages eV — a/S./n, 
with n an integer and a = 2 for standard Cooper pair tun- 
neling and a = 1 for single electron tunneling mediated by 
Majorana fermions. Experimentally, Zhang et al. [118| have 
reported evidences for MAR for junctions that are in the stan- 
dard a — 2 regime (as expected for 3D junctions). Remark- 
ably, in this experiment only steps are observed for even n. 
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Magnetic field (mT) 

FIG. 6: Critical current modulation by magnetic field of Josephson junction with topological insulator interlayers. (a) From 1601 . (b) From 
11171 . (c) From 1 120|. While the Josephson junctions fabricated by Qu et al. 11171 (b) have standard Fraunhofer patterns, the junctions from 
Veldhorst et al. L60] (a) and Williams et al. L120J (c) show deviations. Deviations occur both with respect to the field period (which is smaller 
than expected from the effective junction area) and in the shape. 



Current-phase relationship 

Topological protection causes interfaces to be highly trans- 
parent due to Klein tunneling l!T4l [87)1 and absence of 
backscattering for perpendicular incidence. This protection 
can change the standard current-phase relationship in Joseph- 
son junctions from sin(0/a) with a = 1 to a = 2. The sys- 
tems is then gapless at zero energy, and there exist a Majorana 
bound state when the phase difference across the junction 
equals tt. Many theoretical proposals have been put forward 
to reveal this current-phase relationship. These proposals are 
based on the ac and dc Josephson effects and superconducting 
quantum interference devices with topological interlayers. 

The ac Josephson effect has been observed by Veldhorst 
et al. fSO], see Fig. |5] These Josephson junctions show 
clear steps when irradiated with microwaves, demonstrating 
the Josephson nature of the supercurrent and from the spacing 
between the steps it can be concluded that these measurements 
are in the a = \ regime. Future strategies to enhance the frac- 
tion of sin(0/2) tunneling involve non-equilibrium measure- 
ments 1.1261 and the suppression of non-perpendicular trajec- 
tories in the junctions 0127L 

The dc Josephson effect causes a modulation of the super- 
conducting critical ciment in an applied magnetic field. In the 
limit of infinite width and a homogenous cuiTent density dis- 
tribution the magnetic field dependence of the critical current 
is the Fraunhofer sine function: 



/c($0) = /c(0) 



sin (7r$/a<i>o) 
7r$/a<i>o 



(7) 



Here, $o is the superconducting flux quantum. Fig. |6] shows 
typical critical current modulation patterns observed so far 
The Josephson junctions fabricated by Qu et al. Ml 171 are well 
described by Eq. (j7]i with a = 1, while Veldhorst et al. Il60l 
and Williams et al. II 12011 report deviations. 

Deviations in the field dependence can result from flux fo- 
cusing effects and geometrical inhomogeneities, e.g. pinholes 



result in a slower decay of the side lobes, while lower current 
densities at the edges result in faster decay of the side lobes. 
Furthermore, having a finite width changes the field depen- 
dence, which becomes significant when W ^ L, the situation 
for most experiments so far. The ratio between the length L 
and the width W of the junctions can change both the mag- 
nitude and period of the diffraction pattern 1 128 1. The period 
increases from $o to 2$o for L/W — >■ oo Iil29 l when the 
junction edges are 'open', as in Fig. |7|a). In this scenario, the 
side lobes will decrease more rapidly. 

When specular reflection occurs at the edge of the junction, 
such as in Fig. |7|b), the 2$o crossover occurs at smaller as- 
pect ratiosi 130|: L/W ^ 1. The ratio that is important in this 
scenario is the distance between the Josephson vortex and the 
range of nonlocal electrodynamics, determined by the thermal 
length ^AT. As long as the range of nonlocal electrodynamics 
is smaller than the distance between vortices, the $0 period re- 
mains. For a larger range or strong nonlocality the period be- 
comes 2<i>o due to boundary effects. Including a finite tunnel- 
ing barrier to this geometry mainly causes a sharpening of the 
Fraunhofer pattern and the side lobs are flattened |131 1. This 
flattening is also characteristic for SNS junctions in the diffu- 
sive limit II 1321 . A third situation occurs when the magnetic 
field is parallel to the surface but perpendicular to the current 
as shown in Fig. |7|c). This changes the Fraunhofer pattern 
drastically 11331 . including irregular periods, periods smaller 



than the flux quantum and non-zero minima depending on the 
L/D ratio. The experimentally realized dc SQUIDs by Veld- 
horst et al. 1 11 9 1 based on junctions where the critical current 
field dependency deviates from the standard Fraunhofer pat- 
tern II6OI show standard fluxoid quantization and suggest that 
the junction deviations are due to geometrical effects. 

The observation of Fraunhofer patterns in topological 
Josephson junctions unambiguously shows the development 
of the superconducting proximity effect in the topological in- 
sulator The deviations from the standard Fraunhofer form 
an interesting platform to search for new phenomena in these 



10 




FIG. 7: (a) Geometry of S-TI-S junctions as used in tlie calculation 
of Ref. 1 129] to determine the critical current modulation by applied 
magnetic fields. The dotted arrow shows the path of an electron that 
never goes to the other superconductor but leaves the junction. In 
the middle of the superconductor the angle for which the electron 
still reach the other superconductor is larger than at the edges, (b) 
Geometry as used in the calculation of Ref [ISO]. At the boundaries 
the electrons are specular reflected, (c) The geometry as used in the 
calculation of Ref . 11331 . 

junctions. However, these junctions are nanosized, are often 
ballistic, and multiple bands can contribute to the conduction, 
resulting in complex scenarios. Still, based on these Joseph- 
son junctions SQUIDs can be fabricated that can discriminate 
between geometrical effects and transport mediated by Majo- 
rana fermions. 

dc SQUIDs with unconventional current phase relationships 

Due to the appearance of Majorana fermions single elec- 
tron tunneling occurs and the resulting current-phase relation 
of Josephson junctions can become 47r-periodic. The An- 
periodic current-phase relationship of topological Josephson 
junctions has its signatures in superconducting loops 11341 . 
In the absence of relaxation mechanisms, standard fluxoid 
quantization is doubled. Unfortunately, mechanisms such as 
quasiparticle poisoning and quantum phase slips will drive the 
system easily to standard fluxoid quantization. However, as 
shown by Veldhorst et al. 1 135 1, even in this regime the un- 
usual current phase relation of the individual junctions alter 
the standard SQUID modulation characteristics. This model 
is based on imposing fluxoid quantization on a superconduct- 
ing loop interrupted by two Josephson junctions described 
with the resistivily and capacitively shunted junction (RCSJ) 
model. To further illustrate this model we show the sce- 
nario of full relaxation in Fig. [8] where the junctions have 
been modeled by the BIonder-Tinkham-Klapwijk (BTK) ap- 
proach ||136J . The resulting bound states are described by 

E = ±Y^ ''°''^f+l'^ - Here, Z is the BTK barrier strength 
of the barriers /. These results also represent bound states 
of superconductor - topological insulator Josephson junctions, 
where the factor Z effectively describes the scattering result- 
ing from finite angle incidence with momentum mismatches 
and the appearance of a magnetic gap. In the standard SQUID 
scenario, the Josephson junctions have low transparency (cor- 
responding to Z = 100 in Fig. [8]) and can be described with a 
sin(0) current-phase relationship. In that case the dc SQUID 
has a sinusoidal critical current modulation as a function of 



(a) 




FIG. 8: SQUIDs including transparant junctions, (a) Current-phase 
relationship of junctions with different interface transparancies (Z = 
0, 0.5 and 100). (b) Resulting SQUID characteristics. Dashed lines 
include the peaked current phase relationship as considered in the 
phenomenological model introduced by Williams et al. |T20|, with 
the screening parameter /Bl ~ and 2. These results hold in the 
presence of relaxation mechanisms and are independent on junction 
homogeneities and are thereby strong signatures of unconventional 
current-phase relationships. 

magnetic field and standard fluxoid quantization = h/2e. 
Lowering the barrier strength, increasing the interface trans- 
parency, results in an incomplete critical current modulation 
B135I . This incomplete modulation is a strong signature, since 
it is independent on the junction homogeneity and survives 
up to the regime of energy relaxation (where doubled fluxoid 
quantization is absent |87|). 

We also included the current-phase relationship as imposed 
phenomenologically by Williams et al. [12Q1- This model 
assumes a current-phase relation that peaks when the rel- 
ative phase over the junction equals tt. In the regime of 
small self-induced flux, described by the screening parame- 
ter /3i = the characteristic spikes are also present 
in the dc SQUID characteristics. Increasing smears out the 
spikes, resulting in a triangular critical current modulation, see 
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Fig. |8] We conclude this section by noting that dc SQUIDs 
are ideal candidates to measure the unusual current-phase re- 
lationships of superconductor-topological insulator structures 
even in the presence of junction inhomogeneities and relax- 
ation mechanisms. These dc SQUIDs can be used to test 
whether the critical current field dependences of topological 
junctions are due to intrinsic or extrinsic effects. 
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